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Abstract 

We study bounds on the enstrophy growth for solutions of the viscous Burgers equation 
on the unit circle. Using the variational formulation of Lu and Doering, we prove that the 
maximizer of the enstrophy's rate of change is sharp in the limit of large enstrophy up to 
a numerical constant but does not saturate the Poincare inequality for mean-zero 1-periodic 
functions. Using the dynamical system methods, we give an asymptotic representation of 
the maximizer in the limit of large enstrophy as a viscous shock on the background of a 
linear rarefactive wave. This asymptotic construction is used to prove that a larger growth 
of enstrophy can be achieved when the initial data to the viscous Burgers equation saturates 
the Poincare inequality up to a numerical constant. 

An exact self-similar solution of the Burgers equation is constructed to describe formation 
of a metastable viscous shock on the background of a linear rarefactive wave. When we 
consider the Burgers equation on an infinite line subject to the nonzero (shock-type) boundary 
conditions, we prove that the maximum enstrophy achieved in the time evolution is scaled as 
f , where £ is the large initial enstrophy, whereas the time needed for reaching the maximal 
enstrophy is scaled as log(£). Similar but slower rates are proved on the unit circle. 



1 Introduction 

We consider the initial-value problem for the one-dimensional viscous Burgers equation, 

Ut + 2uux = Uxx, X £ T, t £ 
u|t=o = no, X £ T, 



(1.1) 



where T = [— ^, ^] is the unit circle equipped with the periodic boundary conditions for the real- 
valued function u. Local well-posedness of the initial- value problem (jl.ip holds for uq € H^^j.{T) 
with s > — ^ |8] . The Burgers equation is used as a toy model in the context of a bigger problem 
of how to control existence and regularity of solutions of the three-dimensional Navier-Stokes 
equations [6l [H] . Recent applications of the Burgers equation to the theory of turbulence can be 
found in [I71I2D]. 
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Lu and Doering [T5] considered the question of optimal bounds on the enstrophy growth. The 
enstrophy for the Burgers equation (jl.lh is defined by 

E{u) = ^l^uldx. (1.2) 

Integration by parts for a strong local solution of the Burgers equation ()l.ip in H^^^ (T) yields 



dE{u) 
dt 



f Ux{uxxx - 2uuxx - 2ul)dx = - {ul^ + ul)dx = R{u), (1.3) 



where R{u) is the rate of change of E{u). 

If n S Cpgj.(T), then there is ^ S T such that Ux{i) = 0. Using the elementary bound, 

^( r ('^^^\ , ^ II II ^11 iii/2ii 1,1/2 



and the Young inequality for a, 6 G 



aP e%'i 1 1 

ab<— + , - + - = 1, e>0, 

peP q P q 



the rate of change R{u) in p.3p can be estimated by 



\Riu)\ < -\\Uxx\\l2 + \\Ux\\%^\\Uxx\\]!2 < ^11^x11^2^' = l^'/'W, (1-4) 

provided that p = |, g = 4, and e = \/2. 

In the framework of the Burgers equation (jl.ip . Lu and Doering [15] showed that the bound 
\R{u)\ < CE^/^{u) on the enstrophy growth is sharp in the limit of large enstrophy, up to a 
choice of the numerical constant C > 0. To prove the claim, they considered the maximization 
problem, 

max R(u) subject to E(u) = £, (1-5) 

for a given value of £ > 0. An analytical solution of the constrained maximization problem 
()1.5p was studied in the asymptotic limit of large £ by using Jacobi's elliptic functions. We note 
that the bound (II. 4p is achieved instantaneously in time and it may not hold for solutions of the 
viscous Burgers equation (jl.ip for a finite time interval. 

Ayala and Protas [2] reiterated the same question on the validity of bound (jl.4p integrated 
over a finite time interval. The energy balance equation for the Burgers equation (II. ip is given 
by 

K{u) = - [ u^dx ^^j^^ = f ^(^^^ _ 2uux)dx = -2E{u). (1.6) 

2 Jf dt J J 
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If bound ()1.4p is sharp on the time interval [0, T] for some T > 0, then integration of the enstrophy 
equation (ll.Sp imphes 

E^/\u{T)) - E'l\u,) < \ l^iumt = \ [K{uo) - K{u{T))] . (1.7) 

The Burgers equation (jl.ip maps the set of periodic functions with zero mean to itself. Using 
the Poincare inequality for periodic functions with zero mean, 

K{uo) < -^E{uq), (1.8) 

and neglecting K{u{T)) in (|1.7p . we can obtain 

E{u{T)) < (i^^/3(no) + ^Eino)^. (1.9) 

Note that this bound together with the monotonicity of K{u{t)) < K{uq) implies global well- 
posedness of the initial- value problem (jl.ip for any uq S -ffper 0^) ■ 

Using the extended maximization problem for the global solution of the Burgers equation 
dH]) in Hl,,{T), 

max E(u(T)) subject to E(uo) = £, (1.10) 
«oeHi,,(T) 

Ayala and Protas [2] showed numerically that the integral bound (jl.9p is not sharp even in the 
limit of large £. 

We shall use the notation A = 0{£^) as — )• oo if there are constants C± such that < 
C_ < C+ < oo and C-£^ < A < C+S^. Let T^: be the value of T, where max^^gj^^i^^cj-) E{u{T)) 
is maximal over T € M+. The main claims in [2j are reproduced in Table I. 



Initial condition 


Time 


Enstrophy at T^, 


Energy K at T^, 


A maximizer of (jl.Sp 


o(£:-^-«) 






A maximizer of (jl.lOp 


o(£:-^-^) 







Table I: Enstrophy growth in the Burgers equation from the numerical results in [2]. 

The first line in Table I shows that the instantaneous maximizer of the problem (jl.Sp does 
not saturate the Poincare inequality (jl.Sp and does not lead to large growth of the enstrophy. On 
the other hand, the second line in Table I shows that the bound (jl.9p is not sharp. The bound 
(II. 7p could be sharp if K{uo) — K{u{T^,)) = 0{£^^'^) but the numerical work in [2] reported large 
deviations in numerical approximations of this quantity, 

K{uq) - K{u{T,)) = O(£:0-68±0.25)^ 

which may indicate that the underlying relation may have a logarithmic (or other) correction. 
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In this paper, we shall study further properties of the analytical solution of the constrained 
maximization problem (11.50 . We shall use this solution and its generalizations (see Section 2) as 
an initial condition for the Burgers equation (jl.ip . In particular, we shall address rigorously the 
numerical results of [2]. Our main results are summarized in Table 11. 



Initial condition 


Time T* 


Enstrophy at 


Energy K at 


dlSD and (US]) 


0(£:-2/3iog(f)) 




0(f2/3) 


1^ and (US]) 


0(£:-i/2iogi/2(f)) 


OiS^^ log-3/2(^)) 





Table II: Enstrophy growth in the Burgers equation from our analytical results. 

The analytical results in Table II justify partially the results of numerical approximations in 
Table I. We conjecture that the optimal rate is achieved with 

T, = 0(f-i/2), E{u{n)) = 0{£^/^), K{u{n)) = 0{£), (1.12) 

but we have no proof of this rate at the present time, perhaps, due to technical limitations of 
our method (see Remarks [2] and [3]) . Similarly, we cannot derive an analytical analogue of the 
numerical result (jl.lip and hence, the sharpness of the nonlocal bound (II. 7p remains an open 
question for further studies. 

From a technical point of view, using the dynamical system methods, we prove that the max- 
imizer of the constrained maximization problem (jl.Sp does not saturate the Poincare inequality 
(II. Sp . In the limit of large enstrophy £, this maximizer resembles a viscous shock on the back- 
ground of a linear rarefactive wave. If this maximizer is taken as the initial data to the viscous 
Burgers equation (jl.ip . it does not give the largest change of enstrophy, compared to the case 
when the initial data saturates the Poincare inequality (jl.Sp . On the other hand, if the shock's 
width is used as an independent parameter relative to the background intensity of the linear 
rarefactive wave, the initial data can saturate the Poincare inequality (jl.Sp for large values of £, 
up to a numerical constant, and achieve a faster growth of enstrophy in the time evolution of the 
viscous Burgers equation. 

We note that our construction of the viscous shocks on the background of a linear rarefactive 
wave is similar to the diffusive A^-waves that appear at the intermediate stages of dynamics of 
arbitrary initial data in the Burgers equation over an infinite line [13]. However, these metastable 
states correspond to Gaussian fundamental solutions of the heat equation in self-similar variables 
[31 H], whereas our solutions are obtained on a circle of large but finite period after a scaling 
transformation. Our results still rely on the analysis of the Burgers equation over an infinite line 
subject to the non-zero (shock-type) boundary conditions, where viscous shocks are known to be 
asymptotically stable [TO] flOl [TTl [22] . 

The technique of this paper does not use much of the Cole-Hopf transformation [T] 112] , 
which is well known to reduce the viscous Burgers equation to the linear heat equation. The 
transformation is only used in Sections 5 and 6 to reduce technicalities in the convergence analysis 
for dynamics of viscous shocks in bounded and unbounded domains. At the same time, one 
can think about other techniques to prove the conjecture ()1.12p . which explore the Cole-Hopf 
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transformation in more details. In particular, this transformation shows that solutions of the 
viscous Burgers equation can be studied with the Laplace method for the heat equation. The 
Laplace method is typically used to recover solutions of the inviscid Burgers equations from 
solutions of the viscous Burgers equation in the limit of vanishing viscosity (see, e.g., [211 Chapter 
2] or 116\ Chapter 3] ) . Note that the limit of vanishing viscosity corresponds to the limit of large 
enstrophy in the context of our work. There is a promising way to prove the conjecture ()1.12p 
using the Laplace method after the shock of the corresponding inviscid Burgers equation is formed 
and this task will be the subject of an independent work [9]. 

The paper is organized as follows. Section 2 presents main results. Solutions of the constrained 
maximization problem p.Sp are characterized in Section 3. Self-similar solutions of the Burgers 
equation on the unit circle are considered in Section 4. Section 5 presents analysis of the Burgers 
equation on an infinite line. Evolution of a viscous shock on the background of a linear rarefactive 
wave is studied in Section 6. Proofs of the main results for two different initial data in Table II 
are given in Sections 7 and 8. 

2 Main results 

We shall first reexamine the solution of the constrained maximization problem (jl.Sp . Unlike the 
work of Lu and Doering [15J, we avoid the use of special functions (the Jacobi elliptic functions) 
but use dynamical system techniques to study the limit of large enstrophy £. As a result, we 
obtain the following theorem. Here i7pgj,(T) denotes the restriction of //pgj.(T) to odd functions 
and A = Ol°°{B) as i? — t- oo indicates that ||A||ioo = 0{B). 



Theorem 1 For sufficiently large £, there exists a unique solution G i?pgj.(T) of the con- 
strained maximization problem ^.5\l with n^O) < satisfying 



u^{x) = 4k{2x - tanh(A;x)) + OL°o{k^e~^) 



as A; — )• cxD 



(2.1) 



where k determines the leading order expansions, 



K{u^) 



E{u^) 



R{u^) 



^k^+o{k% 

^k' + 0{k'). 
5 



(2.2) 



(2.3) 



(2.4) 




(2.5) 
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Remark 1 Corollary [I] improves the earlier claims in and in fWj/ based on numerical and 
asymptotic computations, respectively. It shows that the Poincare inequality ^.8\) is not satu- 
rated by the solution of the constrained maximization problem U.5\) . whereas the bound \R{u)\ < 
CE^/^{u) is sharp up to a choice of the numerical constant C > with 

35/3 1 

We shall now consider the time evolution of the Cauchy problem (jl.ip with the initial data 

uoix) = 4:k{2x - f{x)), x£T, (2.6) 
where /c > is a free parameter and / : T — t- M is a fixed function satisfying 

feC\T): /(-x) = -/(x), /Q) =1- (2-7) 

The maximizer of Theorem[T]is represented by (j2.ip . Neglecting the exponentially small terms 
as A: — )• 00, this maximizer can be written in the form (|2.6|) with 



We say that the initial data (j2.6p with (|2.8p represents a shock on the background of a linear 
rarefactive wave, where the width of the shock is inverse proportional to large parameter k. 

The maximizer of Theorem [1] does not saturate the Poincare inequality (jl.Sp in the limit 
/c — 7- 00 (Remark [T]). To allow more flexibility, we can take the initial data (|2.6p in the form, 

_ tanh(;x) 

^^'^^"tanh(//2)' ^ 

where parameter I > may be independent of the large parameter k. Figure [1] shows both 
functions (12. Sp and (12. 9p in the initial data (12. 6p by dashed and solid lines, respectively. If /c = 20 
and 1 = 5, the shock in (j2.9p is much smoother than the shock in 
If uq is given by (|2.6p with ()2.9p . then 



K{uo) = k'^Kil), E{uo) = k'^E{l), (2.10) 



where 



32 8 16cosh(Z/2) 



^' 3sinh(0 
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Figure 1: Initial data ([2:6]) with (fTHI) (dotted) or ([23]) (solid) for A; = 20 and / = 5. 

If I = k, then these expansions yield (j2.2p and (j2.3p up to the error terms. In this case, if 
£ = E{uo) is fixed, then 

k = 0{£^/^), K{uo) = 0{£^/'^), as £ ^ oo. (2.11) 

The function F{1) = K[l)/E{l) is plotted on Figured (left). We can see that there is a 
maximum of the function F at / = /q ~ 3.0, where the maximum is at 

FUq) w 0.025297 < At ~ 0.025330. 

If I is fixed independently of k and \{ £ = E{uq) is fixed, then 

)fc = C(fV2)^ K{uo) = 0{£), as £: ^ oo. (2.12) 

This shows that the initial data (j2.6p with (j2.9p saturates the Poincare inequality (jl.Sp in the 
limit /c — 7- oo up to a numerical constant. Note that the value of the constant prefactor F[Iq) for 
Iq ~ 3.0 is 99.9% of the Poincare constant, compared to the numerical computations in [2j, where 
this prefactor was found from solutions of the extended maximization problem (jl.lOp to be 97% 
of the Poincare constant. 

To apply our method, we shall consider a slow (logarithmic) growth of the parameter / in the 
limit k —7- oo, which yields rates slower than rates (|2.12p . In particular, we shall use the following 
elementary result. 
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Lemma 1 Fix A > and let I := {I + A) log(A:). Then, we have 

k = 0{£^/^log-^/^{£)), K{uo) = 0{£log-^{£)), as £ ^ oo. (2.13) 

Proof. As A;, Z — )■ oo, the leading-order expression for K{uq) and E{uq) are given by 

K{uo) = -k^ + 0^-j, E{uo) = -k'l + 0{e). (2.14) 

With the choice of / := (1 + A)log(/c), we are to solve E{uo) = £, which is equivalent to the 
implicit equation, 

?>£ 

z = x\os{x) + O ix) , z := — — , x:=k'^. (2.15) 

5V ; \ J, 16(1 + A)' ^ ^ 

We look at the asymptotic limit z ^ oo. Setting 

zu 

X :- 



log(z) 

we rewrite the implicit equation (12.15P as the root finding problem, 

(l-M^ + !5li!ll + o(l,,_l = „ (2.16) 

V log(2;) log(z) \log{z)_' 



as z — )• oo. 



We note that /(I) = O (^^^fg^) ^ as z ^ oo and /'(I) = l + (^^^fg^) / 
By the implicit function theorem, there is a unique root in the neighborhood of n = 1 such that 
u = l + p°g(i°gW) ') as z ^ oo. The assertion ([233]) holds by (fm]) and (l2T3|l . □ 



log(2) 

If uo is given by (12. 6j) and (12. 9p . then 



12/ /-^/^ , , _ . 2/2 /-'/s 

16 



/ sech (x)(ix H ^ / sech ix)dx 

Jo ^ ^ tanh3(//2) Jo ^ ' 

sinh^ (x) sech^ {x)dx. 



tanh2(//2) Jo " """" ' tanh''(//2) 7o 

128fc2;3 W/2 



i?(no) = 64/c^ 



tanh^(//2) 

The dependence of R{uo) versus / for different values of k is shown on Figure [2] (right). Since 
R{uo) > for large values of k and /, the enstrophy E{u) grows initially for t > 0. 

We shall construct a solution of the Burgers equation (jl.ip starting with the initial data 
(|2.6p - (|2.7p . We prove that this solution displays dynamics consisting of two phases. In the first 
phase, a metastable viscous shock is formed from the function —4kf{x). In the second phase, a 
rarefactive wave associated with the linear function 8kx decays to zero. We compute the growth 
of enstrophy in two cases: when / = 0{k) and the scaling law ()2.1ip holds and when / = 0{log{k)) 
and the scaling law (I2.13P holds. The following theorem gives the main result of this paper. 
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Figure 2: Function F(l) = K{l)/E{l) versus I (left) and function R{uq) versus I for different 
values of k (right). 

Theorem 2 Consider the initial-value problem with initial data in \2. 6j) and 12. 9\) . Let 

£ = E{uq) be the initial enstrophy. There exists T,, > such that the enstrophy E{u) achieves 
its maximum at u-t = u{-,T^). If I = 0{k) as k ^ oo, then 

T, = 0{£-^/Hog{£)), E{u,) = 0{£), K{u,) = 0{£'^''^), as £: ^ oo, (2.17) 

whereas if I = 0{\og{k)) as k ^ oo, then 

T, = 0{£-^'Hog^'\£)), E{u,) = 0{£'''Hog-^'\£)), K{u,) = 0{£\og-\£)), (2.18) 

as £ ^ CO. 

Remark 2 Several obstacles appear in our method when we consider the case I = 0(1) and 
the scaling law \2.12\l . These obstacles come from the behavior of the solution u{x, t) near the 
boundaries x = ±i as well as from the constraints on the inertial time interval [0, T], during 
which the solution approaches the viscous shock on the background of a linear rarefactive wave. 

3 Proof of Theorem [1] 

To prove Theorem (TJ we obtain a convenient analytical representation of solutions of the con- 
strained maximization problem (jl.Sp . We set v = Ux and look for critical points v E H^^^i^) of 
the functional, 

J{v)= / {vl + v"^ + Xv"^ + ^v) dx, (3.1) 

JT 
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Figure 3: Phase portrait for the Euler-Lagrange equation ()3.4p with £ = 5. 

where A, /i G M are the Lagrange multiphers associated with the following constraints, 

J [ v'^{x)dx = £, [ v{x)dx = 0. (3.2) 
2 J J J J 

The latter constraint ensures that li is a periodic function on T. If v is even in Hp^^{T) and has 
zero mean, then u{x) = v{x')dx' is odd in H'^{T), u(±l) = 0, and hence, u G Hp^^{T). 

The Euler-Lagrange equations associated with the functional J{v) yield the second-order 
differential equation, 

2v"{x) = 3v^{x) + 2Xv{x) + fi, xeT. (3.3) 

Integrating equation ()3.3p over T and using the constraints ()3.2p . we find /j, = —6£. Hence we 
are dealing with the family of integrable second-order equations, 

where / is an integration constant. 

The phase plane of system (13. 4p is given by {v,v') G M^. A typical phase portrait is shown 
on Figure [3l There exist two equilibrium points of the second-order equation (j3.4p . denoted by 
(a_,0) and (a+,0), where 
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I (-X ± \/a2 + ISS) => a_<0<a+. (3.5) 
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Let us define 

I± = -a\ - \a\ + Q£a± I- <I+. (3.6) 

The equilibrium (a_, 0) is a center, whereas the equihbrium (a+, 0) is a saddle point. For / = /+, 
there exists a homoclinic orbit connecting the stable and unstable manifolds of the saddle point 
(a+,0). This orbit can be found analytically, 

v{x) = a+- 4A:2sech2(A;x), A: = ^ ^/ \^ + 18f . (3.7) 

Inside the separatrix loop, there is a family of T-periodic orbits for / G (/_,/_|_), such that T is 
a strictly increasing function of / with 

27r 

— as I I-, and T ^ oo as I ^ I+, (3.8) 

id 



where 



.2 



CO 



-A - 3a_ = \/A2 + 18<? = 4A;2 



If /c > vr, there is a unique Iq £ (/_,/+) such that T = 1 at / = /q- The corresponding 
1-periodic solution v G i7pg^(T) of equation (|3.3|) is a critical point of J{v). Integrating (|3.4|) for 
I = Iq further, we obtain a 1-periodic solution v{x — xq), where xq S T is chosen uniquely from 
the constraints v'{0) = and v{0) < a_ < (which yield an odd u with u'{0) < 0). In this way, 
a critical point of J{v) is obtained and the parameter A needs to be defined by the constraint 
E{u) =£. 

To satisfy the constraint E(u) = £ and to justify the asymptotic expansion (|2.ip . we can use 
the representation, 

v{x) = a+- Ak^yiO, C = kx, (3.9) 
where y(^) is a fc-periodic solution of the second-order equation, 

y" -iy + 6y^ = 0. (3.10) 

The fc-periodic solution of equation (j3.10p . which is called a cnoidal wave, is equivalent to a 
A;-periodic sequence of homoclinic solutions, which are called solitary waves [5j (see also Chapter 
3 in [Ij). This representation uses the theory of Jacobi's elliptic functions. We can obtain an 
equivalent approximation result by using methods of the dynamical system theory pE]. 

Lemma 2 There are > and C > such that for all k G (A;o,oo), the k-periodic solution 
2/(0 of the second-order equation /i3. 1 0\) is close to the solitary wave yoo{C) = sech^(0 with the 
error bound, 

sup \y{C)-y^{0\<Ce-\ (3.11) 
^e[-fc/2,fc/2] 

Proof. We write a /c-periodic sequence of the solitary waves as 

yoiO = ^sech\^ + kn) 
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and decompose the solution y of equation (j3.10p as y{^) = yo(0 + ^(O- After straightforward 
computations, Y satisfies 

LY = F{Y), 

where L = -9| + 4 - 12sech^(,^) and 

F(y) = 6F^ + 12 sech"^ {^ + kn)Y + 6 ^ sech2(^ + A;n)sech2(^ + /c/). 

n6Z\{0} neZ«eZ\{n} 

The operator L : H^{R) L^{R) has a one-dimensional kernel spanned by the odd function 
y'ooiO- The rest of the spectrum of L includes an isolated eigenvalue at —5 and the continuous 
spectrum for [4, oo). Hence the operator L is invertible in the space of even functions. 

Let Tfc := [— and denote the restriction of Hp^^{Tk) to even functions by Hp^j.{Tk). 
Let L : Hp^^{Tj^) — )• Lpgj.(Tfc) be a restriction of L. Then, L is invertible if k is sufficiently 
large, and existence and uniqueness of small solutions Y E Hp^^{Tk) of the fixed-point problem 
Y = L~^F{Y) can be found by using contraction mapping arguments, provided that 



sech^(- + kn) 

nGZ\{0} 



and 



L°°(Tfc) 



Y Y sech^(- + A:n)sech^(- + /c/) 

neZ/eZ\{n} 



converge to zero as /c — )• oo. We check this by explicit computations. There are /cq > and 
Ci) C2 > such that for all k G {ko,oo), we have 



Y2 sech^(- + kn) 

neZ\{0} 



< 



L°o(Tfe) 



neN 



gg-{2n-l)fc 

(1 + e-(2"-i)'=)' 



< Cie- 



(3.12) 



and 



Y Y sech^(- + A;n)sech^(- + kl) 

n&l&\{n} 



< llsech^l 



lL2{Tfc) 



Y sedi^i- + kl) 

iei\{o} 



L°o{Tfe) 



+ 2^ ||sech2(- + A;n)||^^(^^,, 



Y sech^i- + kl) 

iez\{n} 



< 3 ||sech^(-)||^2 



^ sech2(. + 

l&\{0} 



L°o(Tfe) 

+ 4A:i/2 ^ Y ||sech2(. + ^- -r r^^jw^^^j^^ 

neN ieN 

< 2V3Cie-'' + 4Cf /fei/2g-2fc < 026-^ (3.13) 



|sech2(. + A:/)||^^ 
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Existence and uniqueness of Y with the error bound (j3.1ip follow from bounds (|3.12p and (j3.13p . 
the contraction mapping arguments for Y = L^^F{Y) in Hp^^{Tk), and the Sobolev embedding 
of/7p\,(T,)toL- (Tfc). □ 
To define uniquely A in terms of £, we recall the constraints ()3.2p . which are equivalent to 
the scalar equation, 

a+ = Ak / y{£)d^ = 8A;(1 + 0{ke-^)). (3.14) 

J-k/2 

Because a+ is related to k by the exact expression (I3.5p . constraint ()3.14p yields a relationship 
between k and £ given by 

> ?>2k^l + 0{ke-'')] = (3.15) 



A + 4A:2 



9f 



1 + 8^ 



Hence we obtain 



which yields as f — >■ oo, 



and 



' o \ 1/3 

fe=(— +l + C)(£:-i/3) E^oo, (3.16) 



2/3 



a+ = 8/fc(l + 0(A:e-^)) = (48^)^/^0(1). 
It follows from (I3.16P that 

32 

£ = —k^ -2,2k^ + 0{k), as A: ^ oo, 
3 

which recovers (j2.3p . 

Using constraints (13. 2p . the Euler-Lagrange equation (13. 3p . the representation (13. 9p . and the 
constraint ()3.14p . we obtain 

R{u) = 2\£ + ]- I v^{x)dx 

pk/2 

= 2X£ + al{Ak^ - a+) - 32k^ y^i^d^- (3.17) 

J-k/2 



Approximations (i3TT]l . (I3J5I1 . and (l3lB in (ISTTp yield 



256 

i?(n) = —k^ + ©(/c^) as A; ^ oo, 
o 



which recovers 
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To obtain (j2.ip and ()2.2p . we integrate the solution v{x) and write 

u{x) = I v{x')dx' = a+x — 4:kz{kx), (3.18) 



where 



Therefore, we have 



Jo 



kx. 



1 /-i/s /-fc/z 

= al-8ka+ xz{kx)dx + 8k z^{i)di. 
24 Jo J~k/2 



It follows from (j3.1ip that for all large k, there is C > such that z{^) is close to Zoo{C) = tanh(^) 
with the error bound, 

sup \z{i) - z^{i)\ < Cke-\ (3.19) 

Ce[-fc/2,fc/2] 

which recovers (12.ip thanks to the expansion (j3.14p . 

Integrating by parts, we also find the elementary expansion, 

J xtanh{kx)dx = - + J log{l + e~'')dz + 0{k~^e~^) as k^oo. 
Together with (j3.14p and (13.16p . this expansion yields 



K{u) = -k'^ -Wk + 32 log{l + e-')dz + 0{k^e-'')) as A: ^ oo, 
3 Jo 



3 

which recovers (|2.2p . 

To complete the proof of Theorem [H we need to show that the critical point u G H'^^^iT) of 
J{ux) is a maximizer of R{u). However, this follows from the uniqueness of the critical point u*, 
for which R{u^) > 0, and from the fact that R{u) is bounded from above by the bound ()1.4p . 
Therefore, the proof of Theorem [T] is complete. 



4 Burgers equation on the unit circle 

To develop the proof of Theorem [2l we convert the initial- value problem for the Burgers equation 
(jl.ip with initial data (|2.6p - (j2.7p to a convenient form, which separates the decay of the linear 
rarefactive wave and the relative dynamics of a shock on the background of the rarefactive wave. 



Lemma 3 Let uq be given by i2. 6\) - [277\ ). Then, a unique solution u G C(M+, i7j^^,(T)) of the 
Burgers equation is given by 

u{x,t)=p{t){2x-w{^{x,t),T{t))), xeT, teR+, (4.1) 
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where 

and t(;(^,r) is a unique odd solution of the Burgers equation, 

Wr = 2wwi: + Wi:i:, |^|<2(fc-r), rG(0,/c), , . 

w\r=o = f{C/4k), ICI < 2k, ^ 

subject to the boundary conditions w = ±1 at = ziz2{k — r). 

Proof. Although the proof can be constructed by a direct substitution, we will give all 
intermediate details. Recall that the initial- value problem (jl.ip has a unique global solution 
u £ C{R+,H^^^{T)) if the initial data uq satisfies (gSD-dSZ])- Odd solutions in X are preserved 
in the time evolution of the Burgers equation (jl.ip and the Sobolev embedding of H^^^{T) to 
Cper(ir) implies that the boundary conditions u{±^,t) = are preserved for all t > 0. 

Let us look for the exact solution of the Burgers equation (jl.ip in the separable form, 

u{x,t)=p{t) {2x-U{x,q{t))), 

where p{t), q{t), and U{x, q) are new variables. If we choose p = —Ap^ and q = p starting with 
p(0) = po and g(0) = 0, then U{x,q) satisfies the initial- value problem, 

Uq + 4xU^ = 2UU^ + ^U^^, xeT, q>0, 

U\t=o = fix), xeT, ^ ' ' 

subject to the boundary conditions U = ±1 at x = it^. In addition, U is odd in x for any q > 0. 
We find from the differential equations p = —4p^ and q = p that 

1 + 4pot 4 
Solving equation Uq + 4xUx = along the characteristics, we define 

^ = 4x =^ X = Ce^'^, 
dq 

where C is an integration constant. Noting that p{t) = poe~^'^, we define (^,r) by 

oq oq ox 

Integrating these equations, we obtain the substitution, 

Aq _ PQX _ 1 _ Aq^ _ Ppt 



U{x,q) = w{i,T), i=poxe ^ = —— — -, t = -pq{1 - e 



which transforms (|4.4p to (|4.3p . Odd functions U m. x become odd functions w in ^ and the 
boundary conditions C/ = iblatx = ib^ become the boundary conditions it; = ±1 at ^ = 
±\{po- 4t). Setting po = 4A: yields gl]), g21), and ([13]). □ 
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Remark 3 The scaling law U.12\) formally follows from the similarity transformation j^.ip . If 
there exists an inertial range C_ < kt < C+ for some k-independent constants < C_ < C+ < 
OO, where the H^-norm ofw in£^ is k-independent, then in this range, p{t) = 0{k), K{u) = 0{k'^), 
E{u) = 0{k^), where k = 0{8^/^) as £ = E{uq) oo. To prove this claim rigorously, we study a 
convergence of solutions of the rescaled Burgers equation j^.gp starting with the initial condition 
^o(0 = f{i/^k) to the k-independent viscous shock Woa{C) = tanh(^) in a hounded hut large 
domain |^| < 2{k — r) for r = 0{k). For a control of error terms, we have to specify further 
restriction on the initial condition /(^/4/c), which result in weaker statements ^2.11^ and Ii2.18\) 
of Theorem\^ 

The Burgers equation Wr = 2ww^ + w^^ admits the viscous shock, 

t/;oo(0 =tanh(0. (4.5) 

In the initial-value problem, 

{ Wr = 2ww^ + w^^, C e M, r G R+, , . 

\ w\r=o = giC/a), eeK, 

where a > is a parameter, the viscous shock is an asymptotically stable attr actor in the space 
of odd functions g with fast (exponential) decay to ±1 as ^ — )• ±00 |10j . To be able to deal with 
the dynamics of viscous shocks in the initial-value problem (14. 3p on a bounded domain, we shall 
first clarify the dynamics of viscous shocks in the initial-value problem (|4.6p on the infinite line. 



5 Burgers equation on the infinite line 

Let us rewrite the initial- value problem (j4.6p for the Burgers equation on the infinite line by using 
the original (unsealed) variables. 



Ut + 2uUx = Uxx, X G M, t G 
u\t=o = uo, X £ M, 

We impose the nonzero (shock-type) boundary conditions. 



(5.1) 



lim n(x,t) = =fC/oo, t € 1R+, (5.2) 

X— >itoo 

for some Uoo > 0. 

To control the enstrophy on the infinite line, we define 

E{u) = ]- f uldx, R{u) = - f {ul^ + ul)dx, = R{u). (5.3) 

The bound \R{u)\ < ^E^^^iu) on the enstrophy growth is derived similarly to (II. 4p . In the case 

of the infinite line, the maximizer of R{u) at fixed E{u) is not decaying at infinity. On the other 
hand, the result of Theorem [1] becomes now explicit. 
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Lemma 4 The maximization problem, 

max R{u) subject to E{u) = £, (5-4) 

admits a unique odd solution 

u*(x) = — 4/ctanh(/cx), (5-5) 

where k is defined implicitly by £, 

£ = E{u,) = —k\ (5.6) 

and 

Proof. The constrained maximization problem (|5.4p for v = Ux € H^{]R.) yields the func- 
tional, 

J{v) = / [vl + v^ + Xv'^) dx, 

where A G M is the Lagrange multiplier. The Euler-Lagrange equations give 

2v"{x) = 3v'^{x) + 2Xv{x), X G M 

for which the only solution v G is the soliton, 

= -4/c^sech^(/cx), A = 4A;^, 

where k > is arbitrary. Litegrating v{x) with respect to x, we obtain ()5.5p . Integrating v'^, v^, 
and Vx over R, we obtain (|5.6p and (|5.7p . □ 

Remark 4 The Burgers equation i5. 1\) subject to the boundary conditions i5.2\) with Uoo = 4A; 
admits the viscous shock solution, 

Uoo{x) = -4:ktanh{'ikx), k > 0, (5.8) 

which yields R{uoo) = and E(uoo) = 4£'(u=k). 

Remark 5 Using the self-similar variables, 

u{x,t) = -4kw{^,T), i = Akx, T = l6k^t, (5.9) 

the Burgers equation i5. 1\) can be written in the form. 



Wr = 2ww^ + w^^, ^ G M, T G 

W\r=0 =Wo, ^ G M, 



(5.10) 



where uq{x) = —4:kwo{(^). If the boundary conditions i5.2\) are imposed with Uoo = 4A;, then w 
satisfies the boundary conditions lim^^±oow{(,,T) = ±1. 

17 



5.1 Initial condition with / = k 

Let us consider the time evolution of the Burgers equations (jS.ip and (jS.lOp starting with the 
initial condition, 

uo{x) = -4:kt&nh{kx), wo(0 = tanh(^/4), (5.11) 

which is a local maximizer (|5.5|) in Lemma HI Note that the initial- value problem (|5.1U|) with 
initial data (jS.lip is independent of parameter A; > 0. We shall prove the following. 

Lemma 5 For any integer m > 0, there is a constant Cm > such that a unique solution of the 
initial-value problem 115. 10\) with initial data i5.11]) satisfies 



sup 



el«l/2a^(u;(e,r) -tanh(O) < C^e-3"/^ r G M+. (5.12) 



Proof. Using the Cole-Hopf transformation [71 [12] , the Burgers equation (jS.lOp with initial 
data (15. lip admits the exact solution, 

V'(C,r) 

where i/ji^^ r) > is a solution of the heat equation -0,- = ip^^ on the real line M with the initial 
condition ?/^(^, 0) = cosh^(^/4). This exact solution exists in the explicit form, 

^(^,r) = M3 + 4cosh(e/2)e^/^ + cosh(^)e^l , ^eR,TeR+. (5.14) 
8 L J 

As r — )• cxD, the last term in ()5.14p dominates and the solution converges in norm to the 
viscous shock Woo{(,) = tanh(^). To prove this convergence, we rewrite (j5.14p in the form, 

ij{^,T) = ■^e^cosh(^) [l + 4cosh(^/2)sech(^)e-^^/^ + 3sech(^)e-^l . (5.15) 
8 L J 

This representation and the Cole-Hopf transformation (|5.13p yield the compact expression, 

w{t r) = tanh(e) + w{^, r), (5.16) 

where 

~ -3r/4 ^ 2 sinh(e/2) - 4 cosh(e/2) tanh(e) - 3 tanh(Oe-"/^ 
w = e ' sechif — — — — -—7- — . (5.17) 

^' l + 4cosh(^/2)sech(0e-3^/4 + 3sech(e)e-^ ^ ^ 

The bound (j5.12p follows from (j5.17p thanks to the exponential decay in ^ and r. □ 
Corollary 2 Fix S > 0. There exist K > and C > such that for all k > K , we have 

sup|n(a;,t)-Uoo(x)| < ^, for all t > := -^^^tili^^l. (5.I8) 
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Figure 4: Top: The exact solution (|5.16p - (|5.17p for different values of r and in the limit r — t- oo. 
Middle: Enstrophy E versus time r. The dotted curve shows the value of E at the viscous shock 
^oo(C) = tanh(^). Bottom: Enstrophy's rate of^Siange R versus time r. 



Proof. Bound (|5.18p follows from (j5.12p with m = after the transformation (|5.9p . 



□ 



Remark 6 We actually have shown convergence of u to in norm for any s > 0, which 
follows from the exponential decay in ^ and r. 

Figure m shows the exact solution w{S,,t) for different values of r together with the rescaled 
enstrophy E = w'^d^, and its rate of change R = 2 J^{w^ — u)|^)d^ versus r. The relevant 
integrals are approximated by using the MATLAB quad function. The enstrophy E is a mono- 
tonically increasing function to the value Eoo = |, whereas the rate of change R is initially 
increasing and then decreasing to i?oo = 0. Note that by Lemma O we have exponentially rapid 
convergence ||i(;(-,r) — it;oo||j/i — )• as r — oo. 

lf£ = E{uo) = ^k^, then bound ([5T8|) shows that 

E{u^) = 0{£) and = 0{£-^/'^ log{£)) as £: ^ oo. (5.19) 

The asymptotic result ()5.19p is included in the scaling law (12.17P of Theorem [2j We have not 
achieved the maximal increase of the enstrophy with the initial data (jS.lip given by the instanta- 
neous maximizer of Lemma HI To achieve the maximal growth of the enstrophy, we shall modify 
the initial data uq. 

5.2 Initial condition with arbitrary / 

Let us consider a more general initial condition, 

uo(x) = -4A;tanh(/x), wo{£,) = tanh{^ / a) , (5.20) 

where o = ^ is the only parameter of the initial- value problem (jS.lOp with initial data (|5.20p . 
For analysis needed in the proof of Theorem [2l we need convergence of the solution w{£,, r) to 
the viscous shock Woo{C) = tanh(^) in the iJ^-norm as the time r gets large, and the smallness 
of w{^,t) — WooiC) at large values of ^ for all times r > 0. Both objectives are achieved in the 
following result. 

Lemma 6 Fix 6 > 0. There are constants A > and C > such that for all a > A, a unique 
solution of the initial-value problem /i5.10\) with initial data 115. 20]) satisfies 

sup|u;(e,r) -tanh(e)| < 3.,^ for all r > hi + 6falog{a) (5.21) 

Furthermore, for any integer m > 0, there are constants Am > and Cm > 0, such that for all 
> ^m; have 

\d^{w{tT)-tanh{0)\<^^^ for all |^| > 1(1 + 5)^0 log(a) and r > 0. (5.22) 
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Proof. Using the Cole-Hopf transformation (j5.13p . we rewrite the initial- value problem 
(jS.lOp with initial data (I5.20p in the form, 

^\r=Q = cosh'^(^/a), e e M. ^ ' 

The initial data can be represented in the form, 



The function is bounded by its value at ^ = 0, 

1 < 4>a{i) < 2", e e (5.24) 

but the upper bound diverges quickly as a — t- co. On the other hand, for large values of ^, the 
function (j)a decays to 1 exponentially rapidly. 

Fix 5 > and define '■= ^-2^alog(a). For sufficiently large a » 1, there is C > such that 

< Mi) - 1 = e'^i°g(i+^""''^'^) - 1 < Cae-2|«l/'^ < |e| > ^ (5.25) 



a 



The upper bound in (|5.25|) converges to zero as o — )• cxd. 
To bound derivatives of (/>a in ^, we note that 



which implies that 



Writing ([QB]) as 



b'^{^) = -2(1 + e-2?/«)«-ie-2«/^ ^ > 0, (5.26) 



sup |</.',(e)l = 2^ sup |<i6',(OI < (5.27) 



d 2e~^5/" (f 1 /f\ 

di^°^^"^^^ = -TT^' ^ dc^^°^^"(^) = r^^^ UJ' 

and using ()5.24p and ()5.25p . we obtain by induction that, for any integer m > 0, there is Cm > 
such that 

sup \d^MO\ < CrrX, sup \d^MO\ < (5-28) 

The solution of the initial- value problem for the heat equation ()5.23p is obtained explicitly as 

V'(^,t) = —1= r e-^+\^\Uv)dv 



2"^/47rr J-oo 
1 



2"0F 



e"^Va(e + 2t + 2V^z)(iz + e-^ / e"^Va(e - 2r + 2V^z)(iz 

2+ J — CO J 
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where 

'^ = ^- 

We shall rewrite the solution in the equivalent form, 

V'(e,r) = ^e-cosh(e)[l + ^(C,r)], ^(C,r) := ^' (5.29) 

where 

/oo 
e"^Va(e + 2r + 2V^z)dz- V^, (5.30) 

e"^Va(^ - 2r + 2V^z)(iz - V^. (5.31) 

-oo 

Using the Cole-Hopf transformation (|5.13p and the explicit representation (j5.29p . we write the 
solution of the Burgers equation (|5.1U|) in the form, 

«;(C,T)=tanh(C)+*(C,T), ^(e, r) = ^^^|^, (5.32) 

where 

f^ix,(c \ (g, t) + e-^d^jj- (g, r) 1 .2(c\( i (c \ i ic \\ 

h^KM = ^ _ + --^sech {^){iP+{C,r)-ip-{tT)). 

V^(e? + e ^) 2^y^T 

We first prove bound (j5.2ip . Because ■0 is even in ^, we can consider ^ > without the loss 
of generality. To analyze we define 

L _ g-Ca + 2r 
2^- 2V^ ' 

so that z > —Za corresponds to > .^a hi the argument of 4>a{0- For any ^ > and t > Ta ■= 
(1 + 5)^a, we have 

2Ta-Ca ^ 1 + 26 



which shows that — t- cx) as a — t- oo. The term r) can be split into the sum of three terms 

/ e-''dz-V^]+ e-^'(0a(e + 2r + 2V^z)-l)dz 

/ — Za 
e"^' -/>«(? + 2r + 2^/Tz)dz = 1 + 11 + III. (5.33) 

For all ^ > and t > r^, we obtain 



T I -2^ , / 1 -2^ y/2 _(l±2£)!aW(a) 

/ = / e dz < e " < p 8 "^°gw 

2za (1 + 25)Valog(a) 
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and 



/// < 2" / e^'-" dz < e-^- < , §" 

-z+ (l + 2(5)^alog(a) 



where the upper bounds are exponentially small as a — >• oo and the bound ()5.24p has been used. 
Using bound ()5.25p . we obtain 



/oo 
-Za 



c 



Combining these terms and dropping the exponentially small terms, we infer that for any fixed 
5 > and for sufficiently large a ^> 1, there is C > such that 

C 

sup sup |V+(C,r)| < ,^^s+25^ . (5.34) 

To analyze we define 

-e - + 2r 



z„ := z- 



1^ 2y/^ 



so that z < Za corresponds to ^ < —^a in the argument of (l)a{0- For any < ^ < and t > Ta, 
we have 



Za > -2— —2- > — ^alog(a), 

V^a \/2 

which shows that — co as a — )• oo. The term ip-{(,, t) can be split into the sum of three terms 
i^-i^,T) = ( r e-''dz-V^\+ r e-'\4>a{i-2T + 2^z)-l)dz 

\J—oo / J —oo 

rz— 

+ / e-'^^a{i-2T + 2^/^z)dz = I + II + III. (5.35) 

Jza 

Using computations similar to those for '(/'+(^,t), we obtain for any < < and r > Tq, 

/ < 1 . p-4»l°g(") // < ^ in < ^ . p-f alog(a)+alog(2) 

5v^2alog(a) ' ai5+2&^' 5^2a\og{a) 

Again, for any fixed 5 > and for sufficiently large a S> 1, there is C > such that 

C 

sup sup |V^-(g,r)| < 3g_^^^2 - (5-36) 

r>ra 0<5<5a O 

Finally, the exponential factor in front of ip- yields the following simple estimate: 

sup e~2C^_(^^^) < ^(2« _ l)e-2?« < Ce-(l+5)alog{a)+alog{2)^ ^ > (5 
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where the upper bound is exponentially small as a — > oo. 

Because of the symmetry for ^ < 0, we infer from (15 .34^ . ()5.36p . and (I5.37P that for any fixed 
6 > and for sufficiently large a 1, there is constant C > such that 

supsup|^(e,r)| < (5.38) 



This result shows that ^ is small for r > in the denominator of the exact solution (j5.32p . 

We now proceed with similar expressions for d^ip+{^,T) and d^ip-{^,T). Using the represen- 
tation ()5.33p . we find 

3cl = p~^a 

« 2^ ' 



(5.39) 



and 



1 1 

d^n = ^e-''iMCa)-l) + ^ J e'''d,(t>a{i + 2r + 2^z)dz 

= ^ ze-'' {(PaiC + '^r + 2^z)-l)dz, (5.40) 

1 r 2 2 1 1 f^^" 2 

d^III = e-'+MO)-e-'^M^a) +ir^ + 2t + 2^z)dz 

l^T L J 2^T J_^^ 

= / %e-^Va(e + 2T + 2V^z)dz. (5.41) 



2+ 

Using bounds (I5.24p and (j5.25p . we obtain that for any > and r > there is C > such 
that 

C (1+25)2 (J (J (1+25)2 

Q I < _^e-^-F^"iog{a) \Q jj\ < \d,III\ < ^e-^-^"'°s(")+'^'°s(2) 



where in the computations for d^^II we have used the fact that the function -y/re"'^'^/" is mono- 
tonically decreasing for any t > |, whereas Tq » 0(a) as a — >• cxd. 

Combining all bounds and dropping the exponentially small factors, we infer that for any 
fixed J > and sufficiently large a ^ 1, there is C > such that 

sup|5,V^+(e,r)|<-^=-^g^, r>ra. (5.42) 

Similarly, we obtain 

sup |9gV-(^,T)| < ^-^T^, T>Ta (5.43) 

and 

|ae^_(e,r)| < — / |z|e-^ dz, T>0, i>ia, (5.44) 
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which yields 



sup e-25|5.^_(e,r)| < ^e'^?-^ < _^e-(i+5)-i°gW+-i°g(2) 



, T>Ta- 



(5.45) 



Because of the symmetry for < 0, we infer from (I5.42|) . ()5.43p . and (I5.45|) that for any fixed 
> and for sufficiently large o ^ 1, there is constant C > such that 



sup 



e^ggV>+(g,r) + e-gggV^-(g,T) 



< 



c 



ra 



35+2,52 



, T > Ta- 



(5.46) 



We now need to estimate the difference ip-^-{S,,T) — -(/'-(^t). We can write 



For any |^| < ir, we have 



1 



1 

7t 



(r,-;-2T)^ 

e 4r — e 4t 



— r />oo 



-2t 



(i?7 

^77 



+ / ) Mv + 2r)e 4r sinh ( ^ ) dr/ = / + //. 



< 



r2T 
2a roo 



|J| = I 4>a{2T-r])e " ir' sinh [^^]dr] 



('7-|g|) 



e i-^ dr] < 



r-ICI 



e 4t < — — e 16 



Then, we represent 



// = _y_^(^„(, + 2r)-l)e-^si„h(g)<i„ + -^/^ 



e ' 4^ sinh \ ^ ] dr] 



Ila + Ih 



and estimate 



1 



\Ih\ = ^ / sinh ( 1^") dr^<^ 



e dr] < —=e le 



and 



where we have used again the fact that the function -y/re"^"^/" is monotonically decreasing for 
any r > |, whereas Ta ^ 0{a) as a — t- oo. Because / and 11^ are exponentially small in a for all 
T > Ta, there is C > such that 



sup |V'+(C,t) - V^-(g,T)| < ..... < 

iei<^ 



T > Ta. 



(5.47) 
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On the other hand, for any |^| > ^r, there is C > such that 

sech2(e)|V^+(e,r) - ^_(e,r)| < C2"e-^ (5.48) 

which is exponentiahy small in a for all t > Ta- It follows from (|5.47p - (l5.48p that for any fixed 
5 > and for sufficiently large a ^> 1, there is constant C > such that 

supsech2(e)|V'+(e,r) - < -r^i-w, T > Ta. (5.49) 

Representation (j5.32p and bounds (|5.38p . (|5.46p and (|5.49p yield the desired bound (|5.2ip . 

We can now prove bound (|5.22p for m = 0. Using the exponential factors, we have bound 
(15.371) . which is exponentially small as a — t- oo for all ^ > and r > 0. We can rewrite bounds 
(j5.44p and ()5.48p in the equivalent form: 

e-2«|55V'-(e,T)| <C2«e-2?% e>ea, r>0, 

and 

sech2(e)|V+(^,r)-V~(e,r)| <C2'^e-2«^ ^ > Ca, r > 0. 

Therefore, these bounds are also exponentially small as a — t- oo and we only need to show that 
the terms ip+{(^,T) and 9^V+(C)''") remain small for all ^ > (1 + 5)^a and r > 0. We need the 
extra factor (1 + 5) in ^ > (1 + 5)^^ to ensure that Za is bounded from below by 

e - + 2t ^ ^^a + 2r 



^a- „^ ^ o r > V2^= V'^(l + '^)«log(a), 

where we have used the fact that the function ^JT + reaches minimum at tq = ^ > 0. 

Using the splitting of V'+ in (|5.33p . we infer that the estimates for I and 111 produce expo- 
nentially small terms in a, whereas the estimate for 11 produces an algebraically small term in 
a. As a result of analysis similar to the one for ()5.34p . we find that there is C > such that 

sup sup |V;+(e,r)| <Cae-2(i+^)«"/"<^^. (5.50) 

Using the exact expression (j5.39p and the fact that the function -^e"*-'^^""''^^-'^/^^'^-* reaches 
maximum for r > at 

we obtain 

Q^l < ^ p-^(l+<5)alog(a) 

\/alog(a) 

Similarly, using the exact expression (|5.4ip . we obtain 

\di^lll\ < ^ p-^(l+^)alog(a)+alog(2) 

V^alog(a) 
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Using the exact expression (j5.40p rewritten as 
and the bounds (|5.25p and (|5.27p . we obtain 



\d,n\ < ^ _p-<5(l+<5)alog(a) , C 

'4 ' — /TTTTTT^A-i-i /9 



Combining these estimates together, we infer that there is C > such that 

C 

sup sup \d^tp+{^,T)\ < (5.51) 
r>0 5>(i+5)g„ a ^ 

Using (I5.50p and (I5.5ip . we obtain the desired bound (15.22P for m = 0. Bound (I5.22P for 
m G N follows by recursion from the higher-order derivatives of the exact solution (j5.32p . the 
decompositions (|5.33p and (j5.35p . and the bound (|5.28p on the higher-order derivatives of the 
function (j)a- □ 

Corollary 3 Fix 6 > and let a = There exist A > and C > such that for all a > A, 

sup|u(x,t) -Uoo{x)\ < ^3^+^2^1/2 ' for all t>%:= ^ 3^2^^^"^*'"'' ' ^^'^'^^ 



Proof. Bound ()5.52p follows from (j5.2ip after the scaling transformation (j5.9p . □ 

Remark 7 We can also prove convergence of u to Uoo in norm but it requires some decay as 
1^1 —7- 00 and an extension of bound \5.21\) for the first derivative in ^, which we do not establish 
in Lemma\^ Note that bounds 15.21\) and 15. 22\) will be applied on large but finite intervals in S^. 

Figure [5] shows the exact solution r) with a = 10 for different values of r together with the 
rescaled enstrophy E = ur^d^. The integrals in the exact solution (15.32P were approximated 
by using the MATLAB quad function. The behavior of w and E looks similar to the case a = 4 
shown on Figure |4] but it takes longer for E to approach to the limit £^00 = | from the initially 
smaller value Eq = 

Explicit computation with the initial data (15.20p yields E{uq) = ^k'^l, whereas we recall that 
E{uoo) = ^k^- If I = C(l) as A: 00, we have k = 0(£}^'^) and hence 

E{u^) = 0(£^l'^), T^ = 0{e-^l'^\o^(E)), as E ^ 00. (5.53) 

If / = (1 A)log(A;), then Lemma □ implies that k = 0{£}l'^\o^'^^'^{e)). From ([532]) . we 
have 

* 8{l + A)k 



1 + 



log log(fc) 
log(A;) 



(5.54) 
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Figure 5: Left: The exact solution (j5.32p with o = 10 for different values of r and in the limit 
r — oo. Right: Enstrophy E versus time r. The dotted curve shows the value of E at the viscous 
shock Woo{S,) = tanh(^). 



Hence we have 

E{uoo) = 0{£^/^log-^/^{£)) and = 0{£-^^'^\og^/^{£)) as £ ^ oo. (5.55) 

The asymptotic result (I5.55P is included in the scaling law (12.180 of Theorem [2l On the other 
hand, the maximal rate (j5.53p can not be achieved, because if I = 0(1), then a = 0{k) and 
the lower bound on the time r > ^(1 + 5)^alog(a) in bound (j5.2ip is beyond the time interval 
r G (0, k) in Lemma [3] as k — )• oo. 



6 Dynamics of a viscous shock in a bounded domain 

By Lemma [3l the initial- value problem (ll.ip is replaced by the equivalent problem, 

^ Wr = 2ww^ + w^(., |C|<2(/c-t), rG(0,A:), , , 

1 w\r=o = wo, lel < 2k, > 

subject to the initial condition wo{^) = f{^/Ak) and the inhomogeneous boundary conditions 
w = ±1 at ^ = ±.2{k — t), where / is given by either (12.8P or (12. 9p . 
Let us define the norm 

/ .2{fc-r) 

11.^11^2 := / \w{i,r)\^di\ . (6.2) 

\J-2{k-r) J 
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An approximate solution of the initial-value problem (j6.ip can be thought in the form, 



d 

Ws,pp{C,T) = ^logV'app(^,r), 



(6.3) 



where V'app is a solution of the homogeneous heat equation on the real line. For / in either ([2 
or (j2.9p , the solution is constructed in Lemmas [5] or [6l respectively. We shall assume that 



d_ 



< iJ7^app(e, r) < 1, CeR, r G M+, 



(6.4) 



which is a reasonable assumption for a monotonic transition from t«app,o(0 = tanh(^/a) with 
either a = 4ora— T-ooasfc— T-ooto the viscous shock WooiO — tanh(^). 

We write t(^app(C> ''") = tanh(^) + w{£,,t) and assume that for any fixed Cq G (0,1), and for 
sufficiently large A: ^ 1, there is a small > such that 



sup [ \w{2{k — t),t)\ + 

re[0,Cofc] 



d_ 



w{2{k-T),T) 



+ 



w{2{k-T),T) 



< 



k ' 



(6.5) 



where Cfc — t- as /c — )■ oo. 

Furthermore, we assume that for sufficiently large k ^ 1, there is a small > such that 



sup {\woiO - W^pp^oiOl + \ - Wapp,o(0|) < ^ 



(6.6) 



where — t- as A; — )• oo. 

Using the Cole-Hopf transformation, we rewrite the initial-value problem (j6.ip in the form. 



w 



^r = %, |e|<2(fc-r), TG(0,A;), 
^^^=0 = ^0, |C|<2/c, 



(6.7) 



subject to the initial condition -00 = ^'^Y>{Jq'^q{,^)^^^ aiid the Robin boundary conditions 
■0^ = ib^/; at ^ = ±2 (A; — r). Using the decomposition 

V^ = V'app(l + *), (6.8) 

we find the equivalent initial-value problem. 



= ^55 + 2u'app^'g, 
^'|r=0 = *0, 



subject to the initial condition, 

Mi) 



^o(0 



V'app,o(6 



1 = exp 



\i\<2(k-T), TG{0,k), 
lei < 2k, 



[woiC') - w;app,o(^')] d(,') - I 



(6.9) 



(6.10) 
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and the Robin inhomogeneous boundary conditions, 

*g = ±x(r)(l + ^'), C = ±2(fc-T), Te[0,k], (6.11) 

where 

X(r) = 1 - i/;app(2(fe - r), r) = " «)(2(A; - r), r). (6.12) 

By bound (|6.5|) and explicit expression (|6.12|) . for any fixed Cq £ (0, 1) and for sufficiently 
large k ^ 1, there is a small > such that 

sup (|x(r)| + |x'(r)| + |x"(r)|) <^. (6.13) 

re[o,Cofc] 

Because — )• as A; — )■ oo in (j6.13p . the Robin boundary conditions (|6.1ip converge to the 
Neumann conditions as — )• oo. 

By bound (j6.6p and explicit expression (j6.10p . for sufficiently large A: ^ 1, there is a small 
Dk > such that 

sup 1^^0(01 ||^o|lL^<^fc, ll^ollLf +II^oIIl^ (6-14) 

Because Dk — )• as /c — )• oo in (|6.14p . the initial condition (j6.10p converges to as A; — )• (X). 
Using apriori energy estimates, we prove the following result. 

Lemma 7 Assume (fg.^D , /i6.13\) . and ^6A^. Fix Co € (0,1). There exist constants K >1 and 
C > 0, such that for all k > K and r G [0, CqA:], the unique solution of the Burgers equation 
^6. 1\) satisfies 

Ik-w'appllia +\\d^{w-w^^^)\\l, <C{Ck + Dl) (6.15) 

k,T k,T 

where Ck and Dk are defined by i6.13]) and |g. i^[ ). 
Proof. We note the correspondence, 

W = U-app + 

Fix Co G (0, 1) and introduce the energy for the initial-boundary value problem (16. 9p . 

HiT):=hnl2 +h^^\\l. , T€[0,Cok]. (6.17) 
If if(r) is small for all r G [0, CqA;], then there is C > such that 

\\w-w^pp\\l2 <CH{t). (6.18) 

k.T 
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The initial-boundary value problem (j6.9p inherits the local well-posedness of the initial-value 
problem (fTTT]) if ^'o G H^{[-2k,2k]) as long as 

sup |^'(^,r)|<l, r€[0,Cok], (6.19) 

l€l<2(fc-r) 

to ensure that 1 + ^{^, r) > for all |^| < 2{k — r) and r G [0, k]. Recalhng Sobolev's embedding 
of to L°°, we have 

3C>0: sup |^'(C,t)| < C^JH{t), tG[0,CoA:], (6.20) 

ICI<2(fc-r) 

hence the constraint ()6.19p is satisfied if H{t) remains small for all r G [0, C^k]. 

Computations are simplified for the even solutions ^{—^^t) = ^[^,t), which generate odd 
solutions w{—£^,t) = —w{^,t) of the Burgers equation (16. ip . Multiplying the heat equation (16. 9p 
by the solution ^ and integrating in ^ over the time-dependent interval [— ^_|_(r), ^+(r)] with 
^_l-(r) = 2[k — r), we obtain 

^|:ll^lli2_^+2^+W = 2x(r)M/+(T)(l + M/+(r)) + 2u;app(e+(r),r)M/2^(r) 
where ^'+(t) = ^'(^+(r),r). Canceling redundant terms, this expression becomes 

Thanks to property (j6.4p . the rate of change of ||^'||?^2 in r is almost negative definite, except 

k.T 

of the first boundary term, which is small as A; — )• oo. The boundary term is controlled by (j6.13p 
and ()6.20p . We need one more equation to be able to control the energy H{t) for r G [0, CqA;]. 

Differentiating the heat equation ()6.9p in ^, multiplying the resulting equation by ^'g, inte- 
grating in ^ over [— (^+(t), ^+(t)], and performing similar simplifications, we obtain 

= 2x(r)(l + vI/+(r))(%) + (r)-2x3(r)(l + ^+(r))2 

where {^^()j^{t) = ^gg(^+(r), r). For strong solutions of the initial-value problem ()6.9p . we 
obtain by continuity that 



l-m^{T) = {^>^-2^>^) 



= (^'55) + (r)-2x2(r)(l + vI/+(r)). (6.22) 



Hence we have 



+ll<L^dli? , (6-23) 
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The positive last term in (j6.23p is compensated by the negative second term in (j6.2ip in 
the sum of these two expressions. Additionally, we can move the derivative of ^+(t) under the 
derivative sign and obtain 

^H,{r) = 2{x{t) - ^{t))^+{t) - k!{t)^1{t) + 2x='(t)(1 + ^+{T)f 

- ll^ipL^llii. - 11(1 - ^^vv,if'^'^dll^^ (6-24) 

where 

H^{r) := H{t) - x{T){2^+{r) + ^^(r)). 

The last three terms in the right-hand-side of equation (|6.24p are negative thanks to property 
()6.4p . On the other hand, the functions x{'t) ^-nd %'(''") ^-^'^ controlled by ()6.13p . Integrating 
(|6.24p on [0, Cqt] and using Sobolev's inequality (|6.20p and an elementary inequality 2\/ H{t) < 
1 + H{t), we obtain 



Fi(r) - Hi{0) < ^ + H{T'))dT', 



r G [0,CoA:]. 



By the estimate (j6.13p . Sobolev's inequality (|6.2Up again, and Gronwall's inequality, we infer that 
there is a constant C > such that 

H{T)<C{H{0) + Ck), T£[0,Cok]. (6.25) 

By the estimate (j6.14p . we have H{0) < and bound (j6.25p yields control of the first term in 
the bound (|6.15p from the correspondence ()6.18p . 

Let us now control the second term in the bound (I6.15p . Differentiating the heat equation 
()6.9p twice in ^, multiplying the resulting equation by integrating in over [— ^+(r), ^+(r)], 
and performing similar simplifications, we obtain 

= 2(%) + (r)(%5) + (r) - 2xir)i^^^)l{T) + 2{w,^,,^^h{^^)liT) 

where the subscript "+" denotes the boundary value at ^ = ^+(t). For strong solutions of the 
initial- value problem (j6.9p . we obtain by continuity that 

= (%5)+(t) - 2x(t)(%)+(t) + 2Kpp,e)+(^5)+(^) (6-27) 

and 

(^«app,cc)+ = (w^app,r)+ " 2(w;app)+ (u;app,c)+ = 2x(t) (wapp,c)+ " xir). (6.28) 
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Using (felTT) . (I622D, (fOTD . and ([OS]) we can rewrite (fOG]) in the form: 

= -2(2x^(r) + 2x(r)Kpp,e)+-x'(r))(l + ^+(r))^M/+(r) 
-2x\t){1 + M/+(r))2(4x^(r) + 2x{r){w,^^,^)^ - x'{r)) 

To control ||^'^5||^2 , we define 



I 1/2 
I app,< 



+ (2x=^(r) + 2x(r)(u;app,c)+ - x'W) (2^+(r) + ^^(r)) , (6.29) 
where 7 is to be specified below. Thus we obtain 

^H,{t) = 2(vI/+(r) + M/^(r))^(2x='(r) + 2x(r)Kpp,5)+-x'(r)) 

-2x2(r)(l + M/+(r))2(4x^(r) + 2x(r)Kpp,^)+ - x'(r)) 
+8(x(r) - X'(r))^+(r) - 4x'(r)^'^(r) + 8x3(r)(l + ^+{T)f 
+27x(t)M/+(t) + 4x(T)(%)^(r) - H^^H^, - 



-3||(1 - ^app,c)^/^^celli.^^ - 4|k;/J,^^'|li.^^ - 4||(1 - u;app,e)^/2^c||i. 

-^ii<p,€^iiii.-ii(^ + ^app,^««)^/^^ciii^^- 



If we choose 



7 := — inf inf tt'aDD fff > 0, 

re[0,Cofc]Ce[-C+(r),€+(r)] ^^'^^^ 

the last term is negative for any r G [0,CoA;]. By Sobolev's inequality, 

3O0: sup |%(e,r)|2<C7(||%5||2, + ), re[0,CoA;], 

and the smallness of X) the term 

fe,r fc,T 

is also negative. All other integral terms are negative, whereas the boundary terms are controlled 
by the estimate (|6.13p . Sobolev's inequality (|6.20p . and the previous estimate (|6.25p . As a result, 
we obtain 

H2{t) < H2{0) + ^ Hi + Hi{T'))dT' < H2{0) +Ck, T G [0, Cok]. 

k Jo 
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Using ()6.14p . ()6.25p . and (|6.29p . we infer that there is a constant C > such that 



W^idh < C{Ck + Dl), T G [0, Cok]. (6.30) 
Thanks to the exact expression 

and Sobolev's inequahty, 

3O0: sup \'^di,T)\^ <C{\\^^^\\l2 ), TG[0,CoA;], 

bounds (j6.25p and (j6.30p yield control of the second term in the bound (|6.15p . □ 

Remark 8 Recall that Wapp('C,T) = tanh(,^) + w{^,t), where w is small for large values of r. 
// we consider the operator L = 5| + 2tanh(^)5^ on the truncated domain [— CoiCo] subject to 
the Neumann boundary conditions, then the eigenvalues of this boundary-value problem for even 
eigenf unctions are located at the real axis and bounded from above by —1 for any > 0. This 
suggests the asymptotic stability of the zero solution of '^r = but does not imply any good 
bounds on the resolvent operator (Aid — which is needed for the estimates of the remainder 
terms. Moreover, the resolvent operator (Aid — L)^^ may grow exponentially as ^ oo because 
the continuous spectrum of L on the infinite line domain M touches the imaginary axis and the 
zero eigenvalue. See Section 4-4 Scheel & Sandstede U8^ . Apriori energy estimates used in the 
proof of Lemma^ avoid this problem, as well as they incorporate moving boundary conditions in 
analysis of the remainder terms. 



7 Proof of Theorem [2t Case I = 0{k) as k ^ oo. 

We shall consider the initial data (12. 6p and (12. Sp . This corresponds to the choice I = k, which 
represents a more general case / = 0{k) as A; — t- oo. 

By Lemma [3l a solution of the Burgers equation (jl.ip is written in the form 

4:k 16k^t 

u{x,t)=p{t){2x-w{p{t)x,T)), P(0 = Y^Y6^' ^= l + im ' ^^-^^ 
where w solves the rescaled Burgers equation (|6.ip with the initial data 

T«o(0 = Atanh(e/4) ^ V'o(e) = cosh4\^/4), (7.2) 

where 

A = coth(A;/2) = 1 + 
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The approximate solution of the rescaled Burgers equation (|6.ip is given by Lemma [5j It can be 
written in the Cole-Hopf form ()6.3p with 



1 



i^.ppitr) = i (3 + 4cosh(e/2)e^/^ + cosh(C)e^ 



(7.3) 



Note that V'app,o(C) = cosh^(^/4). Assumption (|6.4p is satisfied by the direct computations. 
Assumption (I6.5P fohows from bound (I5.12P of Lemma [5] with 

3C7 > : Ck = Cke^^^^^\ r G [0, CqA;], 

which is exponentially small as A; — t- oo for any Co G (0, !)■ 
The initial condition, 



(7.4) 



^rQ(^) = cosh^(^"^)(^/4) - 1 = exp 



log 



1 - e-^ 

implies that for sufficiently large /c ^ 1, there is C > such that 



cosh 



sup |^o(e)l < Cke-\ ||f o|Il!„ < Ck^'^e-\ + 

|C|<2fc 



fc,0 



Hence assumption (j6.14p holds with Df^ = Ck^^^e ^ for some C > 0. Because is much smaller 
than Ck for all r G [0,CoA;], Lemma [7] yields 



W — W; 



appIlL 



2 +\\di{w-w^pp)\\l2 <Cke-^''-^\ Te[0,Cok]. 



(7.5) 



Applying Lemma O we write an approximate solution of the Burgers equation (jl.ip with 
initial data (12.611 and (12.811 in the form, 



Uapp(x, t) = p{t) {2x - Wi,pp{p{t)x, r)) , p{t) 



Ak 



1 + IQkr 



IQkH 
1 + 16/ct' 



on the time interval 



< mkt < 



Co 



(7.6) 



(7.7) 



1-Co 

Because \p{t)\ = 0{k) as — oo for any t in the time interval (j7.7p . bound (j7.5p and Sobolev 
embedding of to L°° imply that there are constants C > and a > such that 



sup\u{x,t) -u^pp{x,t)\ < Ck sup \w{^,t) -w^pp{^,T)\ < Ck^/'^e "^ 



l€l<2(fc-r) 



(7.8) 



for any t in the time interval (j7.7p . The error bound (|7.8p is exponentially small as A; — )• oo. 
On the other hand, we can consider 



Uoo{x,t) = p{t){2x - Woo{p{t)x)) , p{t) 



Ak 



1 + 16A;t' 



(7.9) 
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where Woo(C) = tanh(^). Bound ()5.18p in Corollary [2] imply that there are 5 > and C > such 
that 

C 



in the inertial range 



sup|'Uapp(x,t) - 'Uoo(a;,t)| < -j, (7.10) 
< 16/tt < -%r, (7.11) 



where 



l-Ceo(*:)- -1-Co' 

3/c 

We note that Coo (A;) — )• as A; — oo. 

In the inertial range (jT.lip . we can compute the leading order approximation of E{u) and 
K{u) from the values of E{uoo) and K{uoo)- Using the representation (|7.9j) . we obtain 

Kiu^) = ]:p\t) + 0{pit)), E{u^) = lp^it) + 0{p\t)), as p(t) ^ oo, (7.12) 

On the other hand, R{u) is approximated from R{uoo) by 

R{uoc) = -8p\t) + 0{p^{t)) as p{t) oo. (7.13) 

Because R{uoo) < 0, the maximum of E{u) occurs at the time t = T^, < iqj^^^^q^^ (fc)^ • Moreover, 
it is clear that E[uoo) and hence E{u) is decreasing for all times t >T^:. 

It remains to prove that = 0{k-'^log{k)) as k — 7- oo or, in other words, that there exists 
Coo (A;) = 0{log(k)/k) < Coo{k) such that T^, occurs inside 



16fc(l-Coo(A;)) " 16fc(l-Coo(A:))' 

If this is the case, then the scaling law (j2.17p of Theorem [2] follows from the error bounds 
and ()7.10p . the triangle inequality, as well as from the previous computations: p{t) = 0{k) as 
k oo and k = 0{E^/^) as £ —f oo. 

To show that T^, = 0{k~'^\og{k)) as A; — >• oo, we compute i?(napp) by using the explicit 
representation (|7.6p with u;app(^,''") = tanh(^) + i(}(^,r), where w is given by (|5.17p . Asymptotic 
computations yield 

R{u,pp) = 2p\t) «r,C-<r,Cc)de-12/(t) / W%,^^di + 0{p\t)) 

Jo Jo 

= 2Np\t)e-^^'^ - Sp\t) + O (e-3"/V(t),/(i)) 

as e"^'^/'^ —7- and p{t) — ?■ oo, where is given by 

N = cosh(^/2) (-28sech3(^) + 139sech5(^) - 120sech'^(0) 

Jo 

/>oo 

+ / sinh(e/2) sinh(0 (26sech^(^) - 30sech6(0) d^. 
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Numerical approximation of the integral shows that A'^ ~ 5.5189 > 0. (The positivity of A'" also 
follows from the fact that R is positive for the approximate solution Wapp shown on the right 
panel of Figured!) Therefore, -R(Mapp) = at the time t = T^, (corresponding to r = r.,, by the 
transformation in (j7.6p ). when piT^,) = 0{e^'^*^^). Since piT^,) = 0{k) as A; — t- oo everywhere in 
(17771) . we have n = C'(log(A;)) < Cq/c or = 0{k~'^ \og{k)) as ^ oo. The proof of Theorem[2] 
for I = k \s now complete. 



8 Proof of Theorem Wt Case / = C(log(/c)) as A: ^ oo. 

We shall now consider the initial data (j2.6|) and (j2.9p . By Lemma [H we fix A > and set 

/ := (l + A)log(A;). 

This choice represents a more general case I = C'(log(A;)) as — t- oo. 

By Lemma [3l a solution of the Burgers equation (jl.ip is written in the form 



•1) 



u{x^t) = p{t) {2x — w{p{t)x,T)) ^ p{t) 



Ak 



IQk^t 



1 + 16/ct' l + 16A;t' 

where w solves the rescaled Burgers equation (j6.ip with the initial data 

woiO = Atanh(e/a) MO = cosh^'^(C/a), 



.3) 



where a = ^ and A = coth(//2). The approximate solution of the rescaled Burgers equation 
(|6.ip is given by Lemma O It can be written in the form (j6.3p with 



V'app = i^e'' cosh(^) 



1 + 



ggV^+(e,r) + e-^V'-(g,r) 
0F(e« + e-€) 



(8.4) 



where ip± are defined by ()5.30p - (l5.3ip . Note that V'app,o(0 = cosh'*(^/a). Assumption (16. 4p is 
satisfied by the monotonicity of the transition from i(^app,o(0 = tanh(^/a) to u)oo(C) = tanh(^). 
Assumption (16. 5p follows from bound ()5.22p with 



k^ 



as long as 



2{k-T) > -(l + 5)2alog(a), r e [0,Cok]. 



(8.5) 



(8.6) 



Constant is algebraically small as /c — )• oo, whereas the constraint (|8.6p yields an upper bound 
on Co, 

1 ...... . {1 + 6? r log(//4)" 



Cn < Cn* := 1 



4k 



;i + 5)^alog(a) = 1 



(1 + A) 



1 



log(A;) 



To ensure that Cq > 0, we require (1 + A) > (1 + 5)^. 
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The initial condition, 

*o(0 = cosh''(^~^)(e/a) - 1 = exp 



2ae- 



log 



l-e-' 

implies that for sufficiently large k ^ 1, there is C > such that 
sup |^o(OI < Cke~' = ^, II^oIIl? „ < - ^ 

l«l<2fc 

C 



cosh 



C 



hence the assumption ()6.14p is satisfied with = i^a^i/2 fo'^ some C > 0. Constant Dk is 
algebraically small as A: — )• oo provided A > ^. Lemma [7] yields 



fc,T 



+ \\d^{w - w, 



fc.T 



< C 



1 ^ 



rG[0,CoA:]. (8. 



Applying Lemma O we write an approximate solution of the Burgers equation (jl.ip with 
initial data (12.611 and (12.911 in the form, 



Ui,pp{x, t) = p{t) {2x - Wi,pp{p{t)x, r)) , p{t) = Y^TiQ^' 
on the time interval 



IQkH 
1 + 16/ct' 



< le/ct < 



Co 



(8.9) 



^.10) 



1-Co 

Because \p{t)\ = 0{k) as /c — t- oo for any t in the time interval (jS.lOp . bound (jS.Sp and Sobolev 
embedding of to imply that there is constant C > such that 

sup |'u,(x,t) - napp(a;,t)| < Ck sup |^«('^, r) - t(;app(^, r)| 

xeT |?|<2(A:-r) 



< C7 



1 ^ log(i+*)/2(A:) 



.11) 



for any t in the time interval ()8.10p . The error bound (j8.1ip is algebraically small as A; — )■ oo 
provided A > | and 6 > 2, instead of our previous constraint A > i. 
On the other hand, we can consider 



Uoo{x,t) = p{t){2x - Woo{p{t)x)) , p{t) 



Ak 



1 + IQkt 



.12) 



where Woo(C) = tanh(^). Bounds (j5.52p and (j5.54p in Corollary [3] imply that there is a positive 
constant C such that 



sup|napp(x,t) - iioo(a;,t)| < 



Ck 



3.5+^2 + 



£log^^+^ 

1/2 logV2(a) ^ ^35+52-1/2 ' 



< 



.13) 
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in the inertial range 



< im < (8.14) 



where 



1 — Coo 1 — Co ' 

" (1 + A) < 

Because of constraint (j8.7p on Co, we reahze that (1 + A) > 3(1 + 5)^ instead of our previous 
constraint (1 + A) > (1 + 6f. 

In the inertial range (j8.14p . the computations of R{uoo), E{uoo) and K(moo) are identical to 
those in Section 6 and yield expressions (f7l2]) and (17131) . In particular, E{uoo) = 0{k^) and 
K( 

Uoo) — 0(A;^) as A; — ^ oo. Modifications of the previous argument show that the maximum 
of E{u) at the time t = T^, occurs for T* = 0{k~^) as A; — )• oo. By Lemma [H we have k = 
0{£^/'^\og'^/^{£)). Substituting this into the previous expressions yields the scaling law (j2.18p 
of Theorem [2] for / = (1 + A) log(/c). The constraints on A and 5 are consistent if 

A>^, 5>2, and (1 + A) > 3(1 + -5)2, 

which can be satisfied, for instance, by the choice A = 48 and 5 = 3. 
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